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Abstract. The experimental and theoretical status of the "near perfect fluid" at RHIC is 
discussed. While the hydrodynamic paradigm for understanding collisions at RHIC is well- 
established, there remain many important open questions to address in order to understand 
its relevance and scope. It is also a crucial issue to understand how the early equilibration is 
achieved, requiring insight into the active degrees of freedom at early times. 



1. Introduction 

In April 2005, the four experiments at RHIC made a joint announcment of the discovery of a 
"perfect liquid" in high energy Au+Au collisions [1, 2, 3, 4, 5]. This discovery has made quite 
an impression in the popular press, resulting in articles in Scientific American (April 2006) and 
Discover magazine (February 2007). It even was called the most important scientific story of 
2005 by the American Institute of Physics [6]. While the press release was somewhat short on 
details, it offered a tantalizing glimpse of the connections with subjects as far afield as black 
holes and string theory, suggesting that they offered possible insight into the strongly-coupled 
system formed at RHIC. 

The foundation of these discoveries was the success of ideal hydrodynamics in modeling the 
experimental data in a robust way. In other words, the system formed at RHIC turns out to 
be best seen as a single system that evolves collectively, rather than an ensemble of individual 
nucleon-nucleon collisions. The goal of this review is to outline the techniques and successes to 
date of hydrodynamics at RHIC, discuss the possible experimental and theoretical limits of this 
interpretation, and show how the current data points to new horizons in both regimes. 



2. The "Perfect Fluid" at RHIC 

In hydro dynamical calculations for RHIC physics, there are typically three basic stages [7, 8]: 
initial conditions, hydrodynamic evolution, and freeze-out. The initial conditions consist of 
defining the energy density in position space, with each fluid cell given an initial velocity (note 
that these are purely classical calculations!). The aspects most interesting for RHIC collisions 
are when the spatial distribution of the initial energy density is highly anisotropic. This is 
because the second stage of the collision takes place while the system remains in local thermal 
equilibrium and evolves according to the evolution equations of hydrodynamics: 



(1) 
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which are closed by the equation of state p = /(e) e.g. p = e/3 for a gas of photons or QCD at 
high temperatures [9]. During this period, gradients in pressure and energy density induce an 
expansion and associated cooling of the system. As the local energy density decreases, so does 
the temperature (T oc e 1 / 4 ). When a fluid element cools to the scale of the Compton wavelength 
of the pion T c h ~ m^, the system is considered to be sufficiently dilute that it decouples or 
"freezes-out" of the evolution and is transformed into hadrons. This occurs as an isotropic 
decay of the available energy density into known hadrons with thermal abundances in the fluid 
rest frame [10, 11]. 

While this physical scenario was invented 
in the early 1950's and studied through 
the 1970's, it was eventually discarded 
as a good model for strong interactions. z 100 

Many rejected the possibility that small 
systems had sufficient time (or, strong enough 
coupling of the different fluid elements) to 
achieve any sort of local thermal or chemical gg 
equilibrium, much less a kinetic equilibirium 
of the final state particles, as suggested 
by Hagedorn in the 1960's [12] based on 

the exponential rise of the hadron mass *i";-7- r -,. i , 

spectrum. Heavy ion collisions are larger and 0.5 1 1.5 2 

generate higher multiplicities, leading to a P T 
general expectation that equilibrium is likely. 

In RHIC collisions, circumstantial hints Figure 1. 
of early-time thermalization shows up in the at ^/snn = 
relative population of various hadron states, 
both as a function of particle mass and 

transverse momentum. The transverse momentum dependence, shown in Fig. 1 from Ref. [13] 
fit by a Boltzmann distribution, has the familiar form of the blackbody spectrum: f(p,m) = 



PHOBOS dN/dpT for charged pions 
: 62.4 GeV. A Boltzmann fit is shown. 



l/(exp[(Vp 2 + m 2 - (i B )/T\ ± 1) 
the place of photons. 



although it is obviously a "strong blackbody" with hadrons in 
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Figure 2. (left) STAR data on hadron 
ratios fit to a thermal model, from Ref. [5] . 
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When the spectrum of each particle species 
is integrated N(m) oc Vj d 3 pf(p,m), it is 
conventional to take ratios of particle yields to 
cancel out the common volume factor. This 
simple "thermal model" [14, 15] depends primarily 
on two parameters, the freeze-out T c h and the 
baryon chemical potential [is (defined such that 
p/p ~ exp(—2fiB/T c h))- The experimental data 
on hadron yields is fit to a model incorporating 
these paramters (and all the known strong decays of 
high mass resonances to the observed particles) to 
estimate T c h and [x b ■ Typical values from recent fits 
are T ch = 177 MeV [5], corresponding to 2 x 10 12 K. 
According to the hydrodynamical scenario outlined 
above, this suggests that all possible hadron states 
are available to the system as it freezes out, so 
the relevant degrees of freedom active before the 
This itself implies that the system 



it may well make sense to assume that 
freezeout were also in thermal and chemical equilibrium 
was substantially hotter than T c ^ for most of the evolution, especially at the very beginning 
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Thus, 10 12 degrees is the coolest the system can be! It should be noted, however, that the 
relevance of a single temperature does not imply that the entire system freezes out in some sort 
of "global" equilibrium. Rather, this temperature seems to reflect the local properties of the 
system induced by the properties of the known hadron spectrum. 

If the system is in local equilibrium just 
before freezeout, then it is certainly plausible 
that it could show hydrodynamic behavior. 
Hydrodynamics implies the relevance of ther- 
modynamics, but of course not the converse 
does not necessarily hold true [10]. At RHIC 
(and other machines colliding nuclei) it has 
been noticed that the event-by-event angu- 
lar distribution is not isotropic in azimuthal 
angle. An "event plane" can be estimated 
from the produced particles, defined as an- 
gle *& R of the long principal axis of the par- 
ticle angles. The azimuthal distribution rel- 
ative to this event plane is found to show a 
strong cos(2[c/) — *f> R ]) dependence, shown in 
Fig. 3, from Ref. [16]. This is especially pro- 
nounced in more peripheral (i.e. lower multi- 
plicity) collisions, where the overlap of the nu- 
clei is shaped like an almond, relative to cen- 
tral (i.e. higher multiplicity) collisions where 

the overlap is essentially isotropic. This leads to a characterization of the event-by-event angular 
distributions in terms of its Fourier coefficients [18]: 




Figure 3. Charged dN/dcj) relative to the event 
plane for for different centralities in ^/snn = 
200 GeV Au+Au collisions from PHOBOS. 
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One of the early and striking results 
from RHIC was the measurement of v 2 - 
conventionally called "elliptic flow" 
function of centrality (e.g. the number of 
participating nucleons N part ), an example of 
which is shown in Fig. 4 from Ref. [19]. The 
data show the qualitative features described 
above, that peripheral collisions have a large 
value of V2, reflecting the large asymmetry in 
the initial distribution of nucleons, while the 
central collisions have a much smaller value. 
Hydrodynamic calculations use estimates 
of various properties (e.g. energy or 
entropy density) of the initial conditions by 
extrapolating from lower-energy data, and 
turn out to compare well to the data for 
Apart > 200 [20]. At lower energies, similar 
calculations often overpredicted v 2 by factors 
of two, so RHIC data is considered striking 
validation of the hydrodynamic approach. 

Comparing the relevant energy and space-time scales implied by the success of the 




50 



1 00 150 200 250 300 350 
< N par,> 



Figure 4. The elliptic flow parameter v 2 vs. 
N par t for y/sjffi = 200 GeV from PHOBOS 
compared with a hydrodynamical model. 



hydro dynamical models, the matter at RHIC is formed under quite extreme conditions. The 
estimates of the formation time relevant for the hydrodynamic calculations were predicted to 
be in the vicinity of To = 0.6 fm/c, or approximately 2 yoctoseconds (10 _24 )[20]. This number 
is far smaller than the time taken a massless particle to traverse the radius of a hadron (r ~ 1 
fm/c) [21]. At this time, the energy density needed to match the data is around e = 30 GeV/fm 3 . 
This should be compared with the energy density of a nucleon in its rest frame, ejv ~ 500 
MeV/fm 3 , which it exceeds by a factor of 60. And in the same sense as T c h is a lower limit 
of the system temperature at early times, these estimates do not preclude even higher energy 
densities at even earlier times. All of this depends on exactly when the system can be said to 
be in local thermal equilibrium, i.e. on the precise value of tq. 

Thus, to summarize, the success of the hydrodynamic models at RHIC suggest that collisions 
there make something that is hotter, denser, smaller, and faster (in the sense of thermalization 
time) than other known liquids. Furthermore, ideal hydrodynamics is inviscid by construction, 
i.e. there are no non-equilibrium processes encoded by the equations, at least until freezeout. In 
this light, the appelation of "perfect fluid" seems quite reasonable. 



3. The Edge of Liquidity: The "Near Perfect Fluid" 
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Figure 5. Frontiers of the hydrodynamic 
paradigm 



However, it does not take much serious inquiry 
to raise a few doubts on this conclusion. No 
viscosity was needed to reproduce the data within 
experimental and theoretical uncertainties. And 
yet, it is not obvious that incorporating viscosity 
to some level would make the models completely 
disagree with the experimental data. Thus, it may 
be a "nearly" perfect liquid instead of a merely 
perfect liquid. But what exactly characterizes a 
"liquid", as opposed to a strongly coupled gas? 
Both are fluids, while liquids are distinguished 
by attractive interactions between the constituents 
that generate a surface tension. Of course some sort 
of attractive interactions generate the bound-state 
hadrons in the final state, but this says nothing 

about earlier times. 

Thus, until the precise form of the interaction can be shown to be uniformly attractive, and 
until it can be demonstrated that the system evinces no viscosity, the appellation "perfect liquid" 
should probably be replaced by the "near perfect fluid" . This does not reduce the importance of 
the RHIC discoveries. Rather it points to future investigations in both experiment and theory 
to elucidate the detailed microscopic properties of this system. In the meantime, we will assume 
the near-perfect fluid interpretation as a paradigm for looking at the RHIC data. The idea will 
be to vary certain parameters in experiment and theory, to see if any major changes take place 
that are qualitatively different from what is observed in the most-central collisions. 

In order to probe the limits of this near-perfect fluid paradigm one can push to the "edge 
of liquidity" (shown schematically in Fig. 5 experimentally (in energy, geometry, or rapidity) 
or theoretically (in thermalization time, length scale, or longitudinal dynamics) to see if the 
qualitative behavior changes. One way to examine the roles energy and geometry simultaneously, 
we can consider the systematics of elliptic flow simultaneously as a function of energy, centrality 
and system size. In the past 15 years, a large data set on V2 has been measured near rj = in 
the center-of-mass frame as a function of energy and centrality [3, 4, 5]. However, it is was only 
realized very recently that the data can be understood to depend on only several variables via 
the construction of "scaling relations" [24, 32]. These hold despite the enormous change in the 



available phase space as the energy increases, as evidenced by the higher multiplicities and the 
wider rapidity distributions. 

Calculations performed with ideal hydrodynamics indicate that the magnitude of v 2 is mainly 
driven by the eccentricity of the initial-state matter distribution [22], where the eccentricity e is 
defined as 



e std ~ 2 2 [A) 

°Y + a x 

is the variance in the direction along the reaction plane and Oy is the variance 
perpendicular to it. 

The top panel of Fig. 6 shows e s td as a function 
of N par t, and illustrates how it goes to 1 for 
small Np ar t and large N par t- While there is no 
principle requiring V2 to be linear with e s td, it 
turns out to be in all numerical calculations. The 
same hydrodynamic calculations suggest that there 
is a maximum value of t>2/e, which has been 
called the "hydrodynamic limit" reflecting that full 
equilibrium is reached [23]. Again, there is no 
known deep principle behind this concept, as it 
arises out of the interplay of transverse pressure 
and longitudinal expansion. However, all of this 
does suggest that v%j€ is an important diagnostic 
variable which removes the "trivial" contribution 
from the overlap geometry and provides access to 
the magnitude of the pressure build-up, and thus 
potentially to the equation of state [8]. 

Experimental data from several collision ener- 
gies, and a range of collision centralities shows 
that the "pressure", as diagnosed via t^/e is a 
simple near- linear function of dN/dy/(S), where 
S = tto~xo~y is the transverse area of the initial mat- 
ter distribution [24] . This suggests that the pressure 
is a function of the entropy density in the trans- 
verse plane (if multiplicity is generally thought to 
be linear with the entropy, as will be discussed more 
later). Of course, it is observed that the data rises 
monotonically with transverse density, or (almost) 
equivalently with N part as shown for Au+Au in the 
lower panels of Fig. 6. It appears that the "hydrody- 
namic limit" has been observed yet, although it has 
been argued by some authors that this indicates the 
approach to equilibration, rather than equilibration 
itself [25]. The current experimental and theoreti- 
cal situation precludes definitive statements either 
way, but at the very least it appears to be universal 
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Figure 6. (top) Standard and participant 
eccentricity for Au+Au and Cu+Cu 
collision vs. N part . (bottom) v 2 , V2/e s u, 
and V2/e par t for Au+Au and Cu+Cu 
collisions at 200 and 62.4 GeV, from 
PHOBOS. 



across changes in both energy and geometry (two 

ways to the edge of liquidity, or fluidity). It also does not seem to be "broken" even in very 
peripheral events or at low energies, where one might expect the low multiplicities to not favor 
local equilibration. 

One way to probe the limits of this picture is to thus go to much smaller systems, Cu+Cu 



collisions for example, which reduce the number of participating nucleons by a factor of three. 
Data on V2 as a function of N par t, shown in the lower panels of Fig. 6, already shows that these 
systems appear very different. The Au+Au data has a maximum at N part ~ 100 and decreases 
to 1-2% at Np ar t ~ 350. Conversely, the Cu+Cu data starts at 4-5% and only decreases to 3% 
in the most central collisions. Already this looks very strange, since by construction, e s td — ► 
for the most central events. How can central Cu+Cu have so much apparent transverse pressure 
when there is presumably no geometric way to generate it? Things are even more confusing when 
plotting V2/e s td as a function of N part for Au+Au and Cu+Cu. There seems to be no connection 
between these two systems. Theoretical calculations typically assume that the matter density 
is a smoothly varying one, described by a Fermi distribution. However, experiments and many 
dynamical models assume that interactions occur by nucleon-nucleon collisions where individual 
nucleons are distributed uniformly throughout the nuclear volume, e.g. as implemented by 
"Glauber Monte Carlo" approaches [26]. The finite number, both in Cu and Au nuclei, lead to 
fluctuations in the matter density which are especially notable in peripheral Au+Au and nearly 
all Cu+Cu collisions. If one assumes the matter density is not defined by an idealized overlap 
zone of the nuclear densities, but by the distribution of the participants themselves, then the 
correct scaling should be achieved using what is called the "participant eccentricity" as first 
defined by PHOBOS [27] 



tpart 



a 



+ 4Ky) 2 



■o 



v& 200 GeV 



^s= 1 30 GeV 



\fi= 62 GeV 




nJs= 20 GeV 



Figure 7. Extended longitudinal 
scaling in Landau's hydrodynamics, 
from Ref.[30] 
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This measures the shape of the region defined by the 
principal axes of the participants, on an event-by-event 
basis, using a Glauber Monte Carlo calculation. This 
quantity has the advantage over e s u of being positive 
definite for all values of impact parameter. Thus, e part 
does not necessarily trend to zero for 6^0, which is now 
understood as being due to fluctuations. 

Scaling V2 for several collision energies, centralities, 
and system sizes by e par t shows a universal scaling of 
V2 Impart with dN/dy/(S), or similarly N par t as shown in 
Fig. 6. This suggests that the geometrical configuration 
of the participants (at least as indicators for the 
spatial location of inelastic interactions) is "frozen in" 
immediately, i.e. it does not require substantial time 
for the deposited energy to locally equilibrate and react 
hydrodynamically. This is of course fully consistent with 
the previous estimates of To, but suggests that there is 
little room for "free streaming" before thermalization [8]. 



4. Initial Conditions 

Rapid equilibration in strong interactions is a physical scenario with a long history [7, 28]. Landau 
and Fermi considered this possibility in the early 1950's when estimating the total multiplicity 
generated in the collision of two nuclei, or even nucleons, by estimating the total entropy 
assuming all of the incoming energy is thermalized in a Lorentz-contracted volume. In the 
hydrodynamic paradigm, the strength of the interactions leaves no natural scales in the problem 
save two, 1) the initial longitudinal size and thus the time scale To oc Az oc I/7 oc and 
the final temperature, or equivalent density, described above [29]. As described in several recent 
references (e.g. Ref. [30], these initial conditions lead to testable phenomenological consequences. 
The multiplicity in the Fermi-Landau approach scales as N oc s 1 / 4 , a power-law behavior that 
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is compatible with a wide variety of total multiplicity data from a variety of systems over 
several orders of magnitude in beam energy. It leads to Gaussian rapidity distributions with 
variance a 2 y = ln(y / i/2m) [29]. Combining these two formulae leads to particle densities that 
approximately scale when viewed in the rest frame of one of the projectiles, as a functionn of 
f] ~ 2/bcam, as shown in Fig. 7. 

Landau's initial conditions and dynamical evolution are in 
some sense very different than the one suggested by Bjorken in 
the early 1980's, where "boost-invariance" was thought to be 
the relevant feature of the dynamics of strong interactions [21]. 
Boost-invariance is in fact a generic feature of hydrodynamic 
solutions, and was noticed by Landau in his hydrodynamic 
calculations [10]. In general, the difference between the two 
is a matter of evolution time and Landau's initial conditions 
have a much different causal structure. Early thermalization in 
the highly compressed initial state leads to strong longitudinal 
expansion controlled by To- Conversely, boost-invariant initial 
conditions lead to no longitudinal pressure gradients, but simply 
v = z/t (a la like Hubble expansion). 

RHIC data do not show any clear signs of boost invariance 
in the final state. Intead they have a strong rapidity 
dependence [31], also seen by the elliptic flow as a function of 
pseudorapidity [32]. Interestingly, they also show an invariance 
when plotted as a function of rj — ybeam, shown in Fig. 8 from 
Ref. [33]. for dN c h/dr], vi and v±, circumstantial evidence 
that particle production in heavy ion collisions undergoes early 
thermalization and thus strong longitudinal expansion. 

It then becomes an interesting question whether all of these 
scaling relations indicate a broader validity of the near-perfect 
fluid paradigm over most the collision evolution. The scaling 
relations hold as a function of energy, centrality, and rapidity. 
They are thus related to the question of thermalization time 
and its effects on longitudinal dynamics. 

This points directly to an even deeper question, that 
of whether or not there is a dynamical length scale in 
the problem: one that reflects the microscopic dynamics 
than simply the boundary conditions. If no such scale 
exists, then there is no natural way to argue that some 
systems, e.g. nucleon-nucleon collisions or e + e~ annihilation 
into hadrons, are "too small" to rapidly thermalize or 
react hydrodynamically. Those systems have been noted to 
have the same multiplicity (i.e. entropy) [34] and freezeout 
temperature [35] as nominally larger systems, and they 
also have been long-known to show extended longitudinal 
scaling [4, 36]. It is conceivable that all strong interactions 
could be understood in part by the near-perfect fluid 
paradigm. The criticism that e + e~ is already understood 
via perturbative QCD calculations can be countered by the 
fact that these calculations seem to show features that are 
parametrically similar to Landau's hydrodynamic model[37, 38]. In the end, all of this will hinge 
on the viscosity of the produced matter, something which will be discussed below. 
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5. Degrees of Freedom 

And yet, even if the near-perfect fluid paradigm is relevant over most of the evolution, one 
major outstanding question remains. What is the fluid made of? How did it come into being as 
a locally equilibrated state of matter? 

Obviously there must be some non- 
equilibrium dynamics to generate the ob- 
served entropy. Nothing in the data discussed 
so far uniquely indentifies which degrees of 
freedom are able to achieve this. The most 
natural assumption would be that the early 
stages are dominated by the dynamics of free 
quarks and gluons, or at least the dynamics of 
quark and gluon fields that are studied using 
Lattice QCD. However, attempts to model 
the existing data on V2 vs. pt for identified 
particles find that a first-order phase transi- 
tion needs to be put in by hand [39]. The 
existing calculations, shown in Fig. 10 do not 
provide sufficient "softening" of the equation 
of state to model the heavier particles which 
are most sensitive to the speed of sound. 
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Figure 10. Data on v% vs. px compared to 
various equations of state, from Ref. [39]. 
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These data have been studied for many dif- 
ferent particles species, which have different 
mass and quark content. Recent PHENIX 
data [40] shows that all of the available data 
on V2 vs. pt lie near one another when plot- 
ted as a function of V2/n q on the Y axis, 
where n q is the number of valence quarks and 
anti-quarks in the hadron, and KE^juq on 
the X axis, where KEt = ttit — m for each 
hadron of mass m. This suggests a scenario 
where freezeout occurs by the recombination 
of "constituent quarks" , particles which have 
a mass of ~ m/n q and the right quantum 
numbers for each hadron. However, the same 
PHENIX paper also suggests that "the scal- 
ing with valence quark number may indicate a 
requirement of a minimum number of objects 
in a localized region of space that contain the 
prerequisite quantum numbers of the hadron to be formed. Whether the scaling further in- 
dicates these degrees of freedom are present at the earliest time is in need of more detailed 
theoretical investigation". If one considers the overall dynamical evolution of the system, with 
many independent stages in principle [38], it is possible that constituent quark scaling is only 
probing the very final stages before formation of the final state hadrons. Thus, it is premature 
to suggest that these RHIC data gives direct information about the early-time formation of a 
truly quark and gluon phase. 

Suggesting that quasi-free quarks and gluons are not the active degrees of freedom at early 
times does not obviously contradict lattice data. In most calculations it is found that e/T , 
which is proportional to the thermodynamically-active number of degrees of freedom, does not 
approach the Stefan-Boltzmann limit even at high temperatures. In fact it falls short of the 



Figure 11. PHENIX data on V2 for different 
particle species, scaled by the number of 
constituent quarks (n q ), from Ref. [40]. 



limit by about 20%. AdS/CFT-based arguments, which model the strongly-coupled Yang-Mills 
plasma as a 10 dimensional black hole, predict precisely a 25% shortfall in the number of degrees 
of freedom [41] . Similar arguments have also shown that the shear viscosity on the gauge theory 
side is proportional to the entropy density, and their ratio has a lower bound [42], i.e. ^ > 4-. 

It is now an active field of investigation to 
understand the precise mechanism which implements 
strong coupling in QCD. Shuryak, Zahed, Brown and 
others have proposed the presence of colored bound 
states in a QCD plasma [43, 44]. This leads to 
larger cross sections and thus stronger coupling in 
the system. And yet, similar lattice calculations, as 
discussed above, find that the fluctuations of quark- 
antiquark fields are not consistent with the active 
degrees of freedom being qq bound states [45, 46]. 

Conventional fluids, shown in Fig. 12 exceed the 
postulated viscosity bound by a factor of ten [42]. 
Figure 12. 4nr]/s for a variety of The success of the hydrodynamic models shown above 
systems, with the postulated lower suggests that the viscosity is almost negligible, possibly 
bound indicated, from Ref. [42]. saturating the bound or even being essentially zero. 

And yet, the precise value of 7]/s has not been 
established using experimental data. It is clear that saturating the viscosity bound would make 
an exciting connection between relativistic heavy ion collisions and a prediction based on physics 
nominally far from the original domain of RHIC physics. 

6. The Future 

Probing the transport properties of the system 
will be the focus of the next generation of RHIC 
experiments. A natural probe of this is the 
transport of heavy quarks. New silicon detectors 
in PHENIX and STAR are being developed to 
measure charmed particles by means of displaced 
decay vertices. By studying the correlation between 
the flow and energy loss of charmed mesons, 
e.g. from recent PHENIX data [47] shown in 
Fig. 13, it is possible to relate their magnitude 
to a charm diffusion coefficient D [48]. This in 
turn allows an estimation of shear viscosity by the 
Einstein relation D ~ 6i]/sT. Already a first 
measurement of this quantity has been performed 
by PHENIX, which suggests the viscosity seen by 
charm quarks falls at most within a factor of 2- Figure 13. PHENIX data [47] on charm 
3 above the viscosity bound. The next generation SU pp re ssion and elliptic flow for Au+Au 
of measurements should substantially improve the collisions, 
experimental precision. 

At the same time, theoretical calculations will also have make equivalent strides. Fully three- 
dimensional hydrodynamical calculations will be required, with full control over initial conditions 
and freezeout (e.g. Ref. [49]). All possible initial conditions, from Landau's to Bjorken's should 
be accessible. This should be coupled with truly systematic studies in order to assign error 
bars to extracted parameters characterizing the initial state, equation of state, and freezeout 
conditions. Then we will have a quantitative handle on just how near we are to the most-perfect 





fluid. 
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